Failure of the a -factor in describing QD laser dynamical instabilities and chaos 
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We show that the long-estabHshed concept of a Unewidth-enhancement factor a to describe carrier-induced 
refractive index changes in semiconductor lasers breaks down in quantum dot (QD) lasers when describing com- 
plex dynamic scenarios, found for example under high-excitation or optical injection. By comparing laser simu- 
lations using a constant a-factor with results from a more complex non-equilibrium model that separately treats 
gain and refractive index dynamics, we examine the conditions under which an approximation of the amplitude- 
phase coupling by an a-factor becomes invalid. The investigations show that while a quasi-equilibrium approach 
for conventional quantum well lasers is valid over a reasonable parameter range, allowing one to introduce an 
a-factor as a constant parameter, the concept is in general not applicable to predict QD laser dynamics due to 
the different timescales of the involved scattering processes. 
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Semiconductor lasers are essential components in optoelec- 
tronics that affect practically every aspect of our daily life. 
There are amble indications that the present technology is ap- 
proaching a stage where fundamental constraints are limiting 
performance. The semiconductor quantum-dot (QD) laser is 
a strong candidate for introducing improvements at the under- 
lying physics level. A very important consideration is dynam- 
ical performance, which impacts almost all applications. 

In semiconductor lasers, the connection between refractive 
index and optical gain plays an important role in determining 
modulation response and dynamical instabilities. Throughout 
literature, this connection is commonly described by assuming 
a linear relationship between changes of the real and imagi- 
nary parts of optical susceptibility with the proportionality 
given by a constant linewidth-enhancement factor a ||T|. For 
example, with optical injection |2| or feedback fSl, laser dy- 
namics depends critically on the phase dynamics of the intra- 
cavity laser field. It is customary to describe this dynamics by 
introducing a constant a into the field equation (UQ. 

For QD lasers, the concept of an a-f actor has been con- 
troversially discussed 1 7 , . Experimental values range from 
near zero L2--12 | to larger values [13 1 up to 60 (HKTSi, with 
different measurement techniques yielding very different re- 
sults 171 [161 [13. Furthermore, frequency chirp under large- 
signal modulation was found to be inaccurately described by 

am. 

Nevertheless, the failure of a in accurately describing QD 
lasers is either not widely known, or ignored. Studies of QD 
laser dynamics have relied on a to describe the amplitude- 
phase correlation |T8^-^22l and experimental measurements of 
the carrier-induced refractive index in QD active media have 
been expressed in terms of the a-factor LI 71 12314251 . In this 
letter, we discuss situations where the use of the a-factor leads 
to incorrect predictions. In particular, for the cases involving 
optical injection or feedback, which are widely studied exper- 
imentally and theoretically by the laser dynamics community 
(SSHIH, we present examples where very different laser dy- 
namics is predicted if an a-factor is used instead of the full 
microscopic description. 

Our model for the QD laser device includes the full time 
dependence of the polarization of the active medium and al- 



lows us to derive the gain and the refractive index in each 
time step without the need to introduce an a-factor. We apply 
a semi-classical approach using Maxwell's equations and the 
semiconductor-Bloch equations |2S|37j. Carriers are injected 
from the bulk into a carrier reservoir (QW), and from there 
into the QDs. A bound ground state (GS) and a twofold de- 
generate (excluding spin degeneracy) excited state (ES) of the 
QDs are considered, both for electrons and holes. The bound 
states are labeled by the index m G {GS^ES}. In order to ac- 
count for the inhomogeneous broadening of the QD transi- 
tions, we distribute the QDs into different subgroups, labeled 
by an index j, with the transition frequency (0^. The distribu- 
tion function f{j) gives the probability to find a specific QD in 
the 7-th subgroup, such that Y^jfU) = 1- We assume a Gaus- 
sian spectral distribution of the QDs with a width (FWHM) of 
60meV. The following dynamic equations describe the time 
evolution of the slowly varying electric field E{t) and the oc- 
cupation probabilities of the QD subgroups p^^, as well as 

the QW and bulk carrier densities per unit area ri^^ and n^^^^, 
respectively (b G {e^h} distinguishes between electrons and 
holes): 
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Here, k is the optical cavity loss rate, K is the injection 
strength of the external optical signal, scaled by the free- 
running electric field amplitude E'^ and the cavity round-trip 
time Tin. The injected signal is assumed to be a monochro- 
matic wave with the frequency (Oinj. The Einstein coefficient 
Wm gives the spontaneous recombination rate in the QDs. The 
charge carrier losses in the QW and in the bulk are given by 
a bimolecular recombination rate and B^^i^, respectively. 
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The sum in Eq.([3]) describes the charge carrier losses in the 
QW due to stimulated recombination between states described 
by an in-plane vector k^^ (the factor of two accounts for spin 
degeneracy). A is the in-plane device area. The pump current 
density is given by /, with the electron charge e^. The com- 
plex amplitude gain g{co,t) is calculated from the adiabati- 
cally eliminated microscopic polarization amplitudes of the 
QD and QW transitions pii and pf^, respectively [i 8l [37]| : 
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where jlm is the respective QD dipole transition moment. T 
is the geometric confinement factor, Ei^g the background per- 
mittivity, /z^^ the height of each QW layer, N^^ the QD 
areal density per QW layer, and the degree of degen- 
eracy of the confined QD states. The transition frequency 
detuning is Aft)^ = {(o}n — (o) and the polarization dephas- 
ing time constant is 72. Here, the QW transitions are as- 
sumed to be detuned far enough from the QD GS transition 
to not influence the gain appreciably, such that only index 
changes due to the QW are taken into account, with g^^ = 
{ia)T)/{ei,g£QAhQ^) [Y^^^Qwpf{t)/E{t)\ . For the polariza- 
tion amplitudes p^ of the QW transitions, an expression sim- 
ilar to Eq. ^ is used, with the corresponding dipole moment 
jlQW and transition frequency detuning of the QW transitions. 
The real and imaginary parts oig{(0^t) correspond to the elec- 
tric field amplitude gain and frequency shift, respectively. 

The collision terms d /dt\^^^ in Eqs. ([2])-(|4]) describe the 
scattering between the different carrier states in the device. 
We account for direct carrier (Auger-) capture processes from 
the QW into the QD GS and ES, as well as relaxation pro- 
cesses between ES and GS, with microscopically calculated 
scattering rates that depend nonlinearly on the QW carrier 
densities |38|. Additionally, we consider carrier-carrier (c- 
c) and carrier-phonon {c-p) scattering between QW and bulk 
states, described within the relaxation rate approximation 
1 39 1, with relaxation rates of 20ps~^ (c-c) and 4ps~^ {c-p) 
(hole rates twice as fast). 

To compare with the QD laser we also model a QW laser 
device by taking only the QW and bulk charge carrier subsys- 
tems into account, which is obtained by eliminating Eq. ^ 
from the dynamic equations. In order to describe nonequilib- 
rium distributions in the QW, we model the ^-resolved carrier 
distribution. The amplitude gain is then given by g^^ . 

In general, the a-f actor is defined as the ratio of the deriva- 
tives of the real and imaginary part of the optical suscepti- 
bility with respect to the total charge carrier number 

N = i:b + + ^nQ" ILj^m fU) VmpQ . Rewritten in 
terms of the amplitude gain g{cOjt) defined in Eq. ([Sj this 
leads to a = - [^Img((o)] / [^Reg((o)] where the relation 
X{co) = 2£iyg/{icor)g{(o) between the optical susceptibility 
and the optical gain g{co) was used. The derivative d/dN in 
the above definition is, however, ill defined. The contribution 
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0.05 ps-i 




48 ps 


T2 100 fs 
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0.60 ^onm 
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IJ.QW 0.50 ^onm 
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0.44 ns"^ 
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0.24ns-i 
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0.15 


h(o 


0.952 eV 


B^iQW) 54ns-^nm2 


Sbg 


14.2 


hQw 


4nm 


10^1 cm-2 



Table 1. Numerical parameters used in the simulation. 



to the optical susceptibility, and thus to a, is different for each 
charge carrier transition in the considered laser system. Near- 
resonant transitions affect mainly the gain, while having little 
effect on the refractive index, whereas off-resonant transitions 
mainly contribute to the index change. To overcome this prob- 
lem, we define an a-factor via the response of the laser to an 
optically injected signal, as also done in several experimental 
studies igOlgll. We define 



a, 
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The above definition overcomes the need to make assump- 
tions about the exact shape of the charge carrier variation dN 
required for evaluating a. Instead, the charge carrier variation 
is determined from the response of the system to the injected 
signal and thus from the intrinsic system dynamics. Other ex- 
perimental setups yield different charge carrier variations and 
thus different variations in g{co) leading to apparently differ- 
ent a. We show below that this issue can be resolved by a full 
dynamic simulation. 

When injecting a monochromatic optical signal from a 
master laser into a semiconductor slave laser, a phenomenon 
called phase locking emerges, where the phase difference be- 
tween the electric field inside the cavity and the injected signal 
becomes constant |42|. The slave laser then emits a constant- 
wave signal with the same frequency as the injected signal. 
Outside the locking range the laser exhibits complex dynam- 
ics, including chaos, excitability and multistability fl3ti45l . 
The parameter space consists of the injection strength K and 
the frequency detuning AVinj = {cOinj — (Ofr)/{27t) between 
the injected master laser signal and the free-running laser fre- 
quency (Ofr of the slave laser. 

The parameter region in the (^, AV;^j)-plane for which 
phase-locking is possible is shown by the colored region in 
Fig. [TJa) and (b) for the QW and QD laser, respectively. It 
forms Arnold tongues. For low injection strengths, the phase- 
locked region is limited by saddle-node bifurcations (black 
solid lines), and becomes limited by Hopf bifurcations (black 
dashed lines) at higher K. Since the injection of the optical 
signal affects the charge carrier distribution of the laser de- 
vice, we expect the a-factor to change when changing the pa- 
rameters K^AVinj due to the change of operating conditions. 
In order to determine an effective a inside the phase-locked 
parameter range, we apply Eq. ^ at each parameter point, by 
slightly increasing K and evaluating the changes in gain and 
index after the transient time. The resulting values for ainj 
are shown by the color code in Fig. [T] It can be seen that the 
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Figure 1. (Color online) Locking tongues in the (^,Av/„j) phase- 
space for (a) the QW laser, (b) the QD laser. The phase-locked re- 
gion is bounded by saddle-node (solid lines) and Hopf bifurcations 
(dashed lines). Black lines correspond to the locking tongue deter- 
mined for the full model and gray lines to a static effective a. The 
color code shows the calculated a-factor inside the phase-locked re- 
gion. / = 2Jfh 



effective a-f actor indeed varies with the injection parameters, 
showing a considerable decrease inside the locking tongue for 
both the QD and the QW laser. To identify the differences 
arising from using a constant a-factor across the parameter 
space, we additionally simulate the QD and QW lasers using 
a gain term corresponding to those used in conventional laser 
models, by assuming a linear relationship between the refrac- 
tive index and the gain, given by g{(o) = (1 — /a) Reg (co) 
where the effective a can be arbitrarily chosen. We use 
a = ainj 3.tK = AVinj = as defined in Eq. ([7]), which should 
approximate the a-factor of the laser under optical injection in 
the best way possible (a = 0.73 in the QD case and a = 1.72 
in the QW case). Surprisingly, the limits of the locking tongue 
shown in Fig. [T] can be very well described by using constant 
a-factors (grey lines) both for the QW and the QD lasers. 

Outside of the phase-locked parameter range, the optically 
injected laser exhibits oscillatory intensity pulsations 1211341. 
These oscillations occur on a timescale comparable to the 
charge carrier scattering lifetimes of QD lasers. We therefore 
expect the dynamics of the gain and refractive index of the 
QD laser to become important here. Fig. [2] shows a compari- 
son of the laser dynamics outside of the phase-locked region 
calculated with the full model (black solid lines) and with a 
constant effective a (red dashed lines). The time- series of 
the intensity in Fig. [2ja) reveal qualitatively similar dynam- 
ics in the case of the QW laser. Also a near-linear relation 
between the real and imaginary susceptibility exists as shown 
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Figure 2. (Color online) Laser dynamics with optical injection out- 
side the locking range, (a) Time-series of the intensity for the QW 
laser with K = 0.3,Av,>,^ = -4.0GHz (left), and for the QD laser 
with K = 0.5,Av/„j = — 4.25GHz (right), with dynamically calcu- 
lated phase dynamics (black solid line) and with constant effective 
a-factor (dashed red line), (b) Trajectory in the complex susceptibil- 
ity plane, shifted by the free-running laser susceptibility Xo for th^ 
QW (left) and QD laser (right). Same parameters as in (a). / = 27^/^ 

in Fig. [2jb), which justifies the assumption of a constant a. 
However, for the QD laser case displayed in the right column 
of Fig. [2ja) the assumption of a constant a-factor leads to 
qualitatively different dynamics (dashed line showing period- 
2 oscillation) if compared to simulations with the full model 
(period-4 oscillations). The reason for this becomes apparent 
in Fig.|2jb), where the trajectory in the complex susceptibility 
plane deviates appreciably from the linear relationship given 
by a = 0.73. Here the independent dynamics of resonant and 
off-resonant states leads to a desynchronization of gain and 
refractive index in the full model. The time evolution of the 
optical susceptibility in this case cannot be described by a sin- 
gle a-factor. Note that the desynchronization of the real and 
imaginary part of the susceptibility is observed for all K in the 
QD laser. 

To shed more light on the differences in the dynamics aris- 
ing from using a constant a, we calculate bifurcation dia- 
grams of the QD and QW laser dynamics outside of the lock- 
ing range, using the full model on the one hand, and a constant 
a on the other hand (see Fig. |3]). By sweeping the injection 
detuning downwards from the phase-locked parameter range 
for a constant K and plotting the intensity extrema of the time- 
series, we can numerically trace bifurcations in the {K^AVinj) 
parameter plane. As discussed before, the assumption of con- 
stant a does not lead to qualitative differences in the dynamics 
of the QW laser (see Fig.jSJa)). The QD laser, instead, reveals 
remarkable differences in the bifurcation structure outside the 
phase-locked region if a constant a is used (see Fig.|3jb)). At 
^ = 0.5 the full model exhibits large region of chaotic dynam- 
ics, followed by inverse period-doubling bifurcations, below 
the saddle-node bifurcation. These bifurcations are missing 
completely if a constant a is used. Also at ^ = 0.7 the con- 
stant a-factor leads to more complex periodic orbits, while 
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Figure 3. (Color online) Bifurcation diagrams of the output intensity 
extrema determined by sweeping the injection detuning Av/„y from 
the phase-locked detuning range towards lower values for (a) the QW 
laser sitK = 0.5 (left), K = 0.6 (right), and (b) the QD laser SLtK = 0.5 
(left), K = 0.7 (right). In (a) and (b) top (red) data are obtained using 
constant a and bottom (black) data using the full model. The saddle- 
node (SN) bifurcation limiting the phase-locked regime is shown by 
the vertical dashed line. J = IJ^h 



the full model predicts only period- 1 oscillations outside the 
locking range. This reveals that the assumption of a linear 
relation between refractive index and gain is not justified in 
QD lasers and will eventually lead to an incorrect prediction 
of the QD laser dynamics. Note that in QW lasers the use of 
a is justified only for sufficiently low injection current and in- 
jection strength (^ ^ 0.7). Otherwise non-equilibrium effects 
become important, invalidating the use of an a-factor also for 
studying complex dynamics in QW lasers. 

So far we showed that the use of an a-factor in QD lasers 
is only valid when discussing steady- states. However, even 
then the a-factor needs to be treated with care as it changes 



for each operation point (i.e., pump current) and experimental 
setup. Apart from the optical injection, common ways of de- 
termining a include the evaluation of frequency modulation 
and amplitude modulation response |46| |47 1 and the evalua- 
tion of amplified spontaneous emission (ASE) spectra |48|. 
From simulations of the different experiments for the QW 
laser, the calculated a is similar for all setups. This is due 
to very fast relaxation processes coupling the resonant (near 
band-edge) and off-resonant (higher k) states. As long as non- 
equilibrium effects can be neglected, i.e., for sufficiently low 
injection strengths, the charge carrier distribution closely fol- 
lows a quasi-Fermi distribution. Therefore a functional re- 
lationship of the ^-resolved carrier distribution on the total 
charge carrier number exists, and thus the derivative d/dN 
is well-defined. In QD lasers, however, we find that differ- 
ent experimental setups yield different values for a, even at 
the same operation point | 37 | due to the considerably slower 
carrier scattering between resonant and off-resonant states. 
Consequently, measurements of the refractive index dynamics 
gathered from one experiment should not be used to predict 
the laser response in a different setup, since the underlying 
charge carrier dynamics may be different. 

To summarize, by applying a semi-classical model to evalu- 
ate the concept of an a-factor in QD and QW lasers, we show 
that the refractive index dynamics in QD lasers is inaccurately 
described by a. We find that in the context of an optical in- 
jection setup it is possible to define an effective a-factor both 
for QD and QW lasers when dealing with cw output, but its 
value varies appreciably with the operating conditions. The 
dynamic response of the QD laser to the injected signal out- 
side of the locking region differs crucially from the dynamics 
predicted by using a constant a, due to the desynchronization 
of gain and refractive index. Thus the concept of a breaks 
down. We expect that result to hold also, e.g., in modulation 
and optical feedback scenarios where the laser emits non-cw 
output. We conclude that the approximations inferred by in- 
troducing an a-factor for the field dynamics are reasonable for 
QW lasers operated close to equilibrium but are too limiting 
for modeling complex dynamic scenarios in QD lasers. 
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